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Doughnut-shaped soap bubbles 
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Soap bubbles are thin liquid films enclosing a fixed volume of air. Since the surface tension is 
typically assumed to be the only responsible for conforming the soap bubble shape, the realized 
bubble surfaces are always minimal area ones. Here, we consider the problem of Ending the axisym- 
metric minimal area surface enclosing a fixed volume V and with a fixed equatorial perimeter L. 
It is well known that the sphere is the solution for V = I and this is indeed the case of a 
free soap bubble, for instance. Surprisingly, we show that for V < ctl? jd'iP', with a « 0.21, such a 
surface cannot be the usual lens-shaped surface formed by the juxtaposition of two spherical caps, 
but rather a toroidal surface. Practically, a doughnut-shaped bubble is known to be ultimately 
unstable and, hence, it will eventually lose its axisymmetry by breaking apart in smaller bubbles. 
Indisputably, however, the topological transition from spherical to toroidal surfaces is mandatory 
here for obtaining the global solution for this axisymmetric isoperimetric problem. Our result sug¬ 
gests that deformed bubbles with V < ctl? jd'iP' cannot be stable and should not exist in foams, for 
instance. 

PACS numbers: 47.55.D-, 47.55.db, 47.55.df 


I. INTRODUCTION 

Soap bubbles have been attracting the attention of 
physics and mathematicians for more than two centuries 
[I]. A soap bubble is a thin liquid film enclosing a given 
volume of air. Surface tension is usually assumed to be 
the only responsible for conforming the bubble surface 
shape, and hence the realized surfaces are always minimal 
area ones. It is well known that the sphere is the solution 
for one of the most celebrated isoperimetric problems: to 
find the minimal area surface enclosing a fixed and given 
volume. Free soap bubbles are known to be spheres. 

We consider here the problem of finding the axisym¬ 
metric minimal area surface with two simultaneous con¬ 
straints: a fixed enclosed volume V and a fixed equatorial 
perimeter L. Since a sphere of radius a is the minimal 
area surface enclosing a volume V = 47ra^/3, it will be 
also the solution for our problem for this volume and 
equatorial perimeter L = 2na. We are mainly interested 
in the cases with L = 2'Ka and V < 47ra^/3, for which 
the solutions may have the shape of a “lens” formed by 
the juxtaposition of two spherical caps of height h < a, 
see Fig. Qa)- The enclosed volume by these lens-shaped 
surfaces are 


U = — {3a^ + h^) , (1) 

whereas their surface area is given by 

A = 27r (a^ -I- /i^) . (2) 

As we can see, for a fixed equatorial perimeter L = 2'Ka, 
one can effectively have arbitrarily small enclosed vol¬ 
umes V by choosing arbitrarily small cap heights h since 
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0 < U < IQk'^ . On the other hand, the surface area A 
will be always bounded from below by a positive value, 
A > LP' 12k. It is clear that for small cap height h, the 
lens-shaped surface cannot be an efficient surface for en¬ 
closing a small volume V with a fixed equatorial perime¬ 
ter L. Surprisingly, in order to obtain the global mini¬ 
mum for this axisymetric isoperimetric problem, a topo¬ 
logical transition is mandatory: from the spherical lens¬ 
shaped to toroidal surfaces. As we will show, lens-shaped 
surfaces of equatorial perimeter L are global solutions for 
our minimal area problem only for volumes V such that 


^ ^ - 6 ^’ 


( 3 ) 


with a ~ 0.21. For V < aL^/dK"^, the axisymmetric min¬ 
imal area surface enclosing a volume V will be necessarily 
doughnut-shaped as that one depict in Fig. [^b), as we 
will see by considering all solutions of our isoperimetric 
problem in the following section. 


II. THE ISOPERIMETRIC VARIATIONAL 
PROBLEM 


Strictly speaking, the isoperimetric problem, dating 
from the antiquity, concerns finding the plane figure of 
maximal area with a given perimeter. In a broader sense, 
however, it includes, for instance, the problem of finding 
the function f(x,y) which minimize a given functional 

s[f] = jj^ C{x, y, /, /a;, fy) dxdy, (4) 

but subjected to integral constraints of the type 



C{x, y, /, fx, fy) dxdy = constant, 


( 5 ) 


where I? is a region of the plane {x, y) and the indices 
X and y denote the respective partial derivatives. Every 



following expression for the area functional 
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a~0.2^ 

(a) 



0 = 60 . 



0 = 100 . 


FIG. 1: Top: lens-shaped surface of minimal area with 
perimeter L and volume V = aL^lGn^, with a ~ 0.21. 
No stable lens-shaped surface with V < should 

exist. Bottom: doughnut-shaped surface of minimal area 
with perimeter L and V = a.1? . Axisymmetric mini¬ 

mal area surfaces with V < aL^fOTv^ are necessarily of this 
type. No doughnut-shaped minimal area surfaces exist with 
V > aL^fGTT^. In both cases, the angle 6 is the internal angle 
of the surface at the equatorial perimeter. 


function here is assumed to be real and smooth. For our 
purposes, let us consider the surface (x, y, z) defined by 
the function f{x,y) = ±z, which we will assume to be 
non-negative and such that f{x,y) = 0 for (x,y) € dV. 

The standard treatment for the isoperimetric problems 
involves the associated functional defined as 

S*[f] = JJ dxdy, C* =£ + AC, (6) 

where A is a constant (the Lagrange multiplier). The 
function /(x, y) that extremizes Q subjected to the con¬ 
straint ([^ also extremizes the free functional Q, i.e. 
f{x,y) is a solution of the Euler-Lagrange equations for 
the associated functional The constant A is to be de¬ 
termined, among all the integration constants, from the 
boundary conditions and the integral constraint ([^. For 
the problem of the minimal area axisymmetric surface, 
one can introduce appropriate polar coordinates (p, 0) 
centred at the surface such that / = f{p), leading to the 


S[f] =27 t [ 
Jv 


while the fixed volume constraint will read simply 


V[f] = 27r / f pdp = constant. 

Jv 


( 7 ) 


( 8 ) 


Clearly, the equatorial perimeter will be given by the 
length of dV. In these coordinates, the Euler-Lagrange 
equation for the associated functional (§ is 



= A, 


which be easily integrated and leads to 


Pf 

/ , — oP +Gi. 

\/l + 2 


( 9 ) 


( 10 ) 


We have two qualitative distinct cases according to the 
value of Cl. For Ci = 0, we have /'(O) = 0, which is 
indeed a regularity condition for axisymmetric surfaces. 
However, and more importantly, in this case there is no 
restriction for the values of p and, consequently, I? is a 
circle. It is quite simple to verify that the solutions of 
(101 for this case are the arcs given by 


{fip)+br + p^ = rl 


( 11 ) 


with rg = 2/A and rg > b > 0. These solutions cor¬ 
responds to the usual spherical cap with basis radius 
= Tq — b^ and height h = tq — b. These caps form 
the lens-shaped solutions for our isoperimetric problem. 

Nevertheless, we have also the solutions with Ci ^ 0. 
Solving (10) for /' and considering the convenient signs, 
we have 


fip) = 


d — x^ 


y(x2-xCj(x^ax-2;2)’ 

where x = Ap, d = 2ACi > 0, and 

^min — a/I d 1, Xmax — 's/1 d 1. 


( 12 ) 


(13) 


The first observation here is the most important one: for 
Cl ^ 0, there will be necessarily restrictions for p, we 
have indeed Xmin < Ap < Xmax- The region T) is not any¬ 
more a circle, but effectively a ring domain. The integral 
(121 can be solved analytically by using elliptic functions. 


but for our purposes we opt to solve it numerically. All 
pertinent details are presented in the Appendix. The sec¬ 
ond additive integration constant will be chosen in order 
to have /(pmin) = 0, with Xmin = Apmin- Notice that 
f'iPmin) diverges, assuring in this way that the juxtapo¬ 
sition of the superior and inferior parts of our doughnut¬ 
shaped surface will be indeed smooth along the interior 
radius. An example of solution /(p) is depicted in Fig. 
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FIG. 2: Solid red line: solution for (12 \ with boundary con¬ 
dition /(pmin) = 0 for d ~ 1.6235, which corresponds to 
Apmin ~ 0.6197, Ap* « 2.2491, and Apmax ~ 2.6197. The 
doughnut-shaped solution depicted in j^b) is obtained by 
the revolution around the vertical axis of the closed curve 
formed by the solution and its reflection on the horizontal 
axis (dashed blue line) for pmin < p < p*. This particular 
value of d corresponds to the toroidal solution with maximal 
enclosed volume. The doughnut-shaped surface is regular ev¬ 
erywhere except on the equatorial perimeter p = p*. 


FIG. 3: Area x Volume diagram for axisymmetric minimal 
area surface with fixed equatorial perimeter L = 2n. The 
solid red line corresponds to the doughnut-shaped solutions, 
with the arrows indicating the direction of increasing d. The 
dashed green line corresponds to the lens-shaped solutions. In 
the detail, the region corresponding to the topological tran¬ 
sition. The maximum volume for the doughnut-shaped solu¬ 
tion is V « 0.869, corresponding to the case depicted in Fig. 

An animation illustrating the transition from spherical to 
toroidal surfaces is available at [2]. 


which corresponds to he doughnut-shaped surface in 
Fig. [3b). The equatorial perimeter of the surface will 
be given by L = 27rp*, where the radius p* > pmin is such 
that /(p*) = 0. 

The constant A can be effectively absorbed by a global 
rescaling. For each value of d > 0, we have a toroidal 
surface. From (13), we see that small values of d corre¬ 
spond to the cases with Smin ~ 0 and ccmax ~ 2. These 
solutions can enclose arbitrarily small volumes, but their 
area is bounded from below by the area of a lens-shaped 
solution with same equatorial perimeter L. On the other 
hand, the solutions with large d, which corresponds to 
large a^min and a:i„axj can enclose arbitrarily small vol¬ 
umes with arbitrarily small surface areas. The situation 
is depicted in the Area x Volume diagram of Fig. [^ 
The solid red line corresponds to the doughnut-shaped 
solutions, while lens-shaped ones correspond to the the 
dashed green line. The doughnut-shaped solution with 
maximal volume corresponds to the case with d « 1.6235 
(depicted in Fig. Qb) and Fig. For any other value of 
d, there are always two minimal area surfaces: one corre¬ 
sponding to small ccmin/a^max (small d), and the other to 
small (xmax — a:min)/a;max (large d). The second one will 
be the the global minimum of the problem, see Fig. We 
see from the diagram that the lens-shaped surfaces are 
effectively the only minimal area solution for our prob¬ 
lem provided that the condition (|3 holds, with a « 0.21, 
which corresponds namely to the minimal area doughnut¬ 
shaped surface of maximum volume. For V < aL^/67r^, 
we see from the diagram that three minimal area surfaces 


coexist, but clearly the global minimum corresponds to 
the case of doughnut-shaped surfaces with d > 1.6235. 
An animation illustrating the topological transition from 
spherical to toroidal surfaces according to the value of the 
ratio V/L^ is available in the Supplementary Material. 

It is interesting to relate the topological transition of 
the minimal area surfaces to the dihedral angle 9 be¬ 
tween the tangent planes at the equatorial perimeter. 
Fig. a depicts the dependence of 9 on the volume V 
for a fixed equatorial perimeter L = 2tt. For the lens¬ 
shaped surfaces, the minimal volume V = aL^/dir^ case 
corresponds to 9 close to 60 degrees. This is the case 
depicted in Fig. Qa). The doughnut-shapes surface with 
the same volume has a larger dihedral angle, close to 100 
degrees (Fig. gb) and Fig. [|). Notice that 

a = 8-^V3 = 0.20577... (14) 

corresponds to the volume of a lens-shaped surface with 
dihedral angle 0 = 60 degrees at the equatorial external 
perimeter. 

Notice that all solutions to our isoperimetrical prob¬ 
lem are indeed constant mean curvatures. This can be 
checked by recalling the first and second, respectively, 
fundamental forms for our surface of revolution gener¬ 
ated by f{p)'. E = F = 0, G = 1 -I- and 
L = -pf'/^/l + r\ M = 0,N = -/7\/l + /'"; and 
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FIG. 4: Dihedral angle 6 between the tangent planes at 
the equatorial perimeter as a function of the enclosed vol¬ 
ume V for axisymmetric minimal area surfaces with equato¬ 
rial perimeter L = 27r. The solid red line corresponds to the 
doughnut-shaped solutions, with the arrows indicating the di¬ 
rection of increasing d, while the lens-shaped solutions are the 
dashed green line. Notice that for d —>■ oo, V —>■ 0 and the 
dihedral angle tends to 180 degrees. The doughnut-shaped 
surface in this limit tends to the usual torus of circular sec¬ 
tion. 


that 


1 / pf \ 

pdpy + /'2 j 


L N 
S + G " 


(15) 


where H stands for the mean curvature of our surface. 
Thus, the Euler-Lagrange equation § is equivalent to 
the constraint of constant H. Since H is constant, we 
can easily evaluate it by taking the point p such that 
f'{p) = 0, where L — Q and N = —f''{p), leading simply 
to 


their solutions accordingly to the ratio V/L^. The typi¬ 
cal lens-shaped surfaces formed by two spherical caps are 
not the global minimum for small enclosed volumes V. 
The global minimum for the axisymmetric case enclos¬ 
ing small volumes corresponds to toroidal surfaces. This 
situation considered here resemble in many ways the clas¬ 
sical Goldschmidt discontinuous minimal area surface of 
revolution limited by two coaxial rings separated by a 
distance £ DE). In our case, for V below a critical value, 
we have two toroidal minimal area surfaces enclosing a 
given volume. One of them has area smaller than the 
corresponding lens-shaped surface, while the other has a 
greater superficial area. In the Goldschmidt case, for £ 
below a certain critical value, we have always two mini¬ 
mal surfaces (catenoids), but only one of them can effec¬ 
tively have a total superficial area smaller than the dis¬ 
continuous Goldschmidt solution. In fact, the diagram 
Area x £ for the catenoids is very similar to our Area x 
Volume curve in Fig. 

Finally, we cannot ignore that a doughnut-shaped min¬ 
imal area surface will be ultimately unstable due to phe¬ 
nomena like shrinking Uli] and Plateau-Rayleigh PIT] 
instabilities. This means that, for instance, if one de¬ 
forms a bubble axisymmetricaly in such a way that its 
equatorial perimeter L is enlarged while its volume V is 
kept constant, the bubble will be inexorably destroyed 
when > Gn'^V/a. If a doughnut-shaped bubble is 
formed, it will probably break apart in smaller daugh¬ 
ter bubbles (see, for similar behavior in another context, 
0). It is not difficult to envisage a non-axisymmetric 
surface with fixed external perimeter L, enclosing a vol¬ 
ume V < aL^/6TT^ and with area A smaller than the 
area of our toroidal surface. Gonsider, for instance, a 
sphere with a handle as shown in Fig. The external 


H = -lrip), ( 16 ) 

which is positive for our toroidal surfaces. The mean cur¬ 
vature can be expressed also as 2H = , where 

i?i and i ?2 are the radii corresponding to the principal 
curvatures. The d —>■ oo (V —> 0) limit of figure @ , for 
which the external dihedral angle tends to 180 degrees, 
is an usual torus of circular section for which Ri R 2 , 
assuring in this way that 2H « R^^ = constant. 



III. FINAL REMARKS 

We have shown that the isoperimetric problem of find¬ 
ing an asixymmetric minimal area surface enclosing a 
fixed volume V and with a fixed equatorial perimeter 
L exhibits a rather unexpected topological transition in 


FIG. 5: Sphere with a toroidal handle: by shrinking the 
handle smaller radius, we have a surface with fixed exter¬ 
nal perimeter L, enclosing a volume V < aL^IG- k^, and with 
area A smaller than our axisymmetric toroidal surface. Ax¬ 
isymmetric surfaces are not global solutions for the problem. 
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handle guarantees the constant perimeter constraint. By 
shrinking its smaller radius, its contribution for the total 
V and A will be arbitrarily small, and hence they will cor¬ 
respond to the sphere, which is known to be the global 
solution for the problem and will certainly encapsulate 
a given volume V with surface area A smaller than our 
doughnut-shaped surface. This kind of non-axisymmetric 
surface might arise from a Plateau-Rayleigh instability, 
where perturbations with wavelength greater than the 
smaller dimension of the toroidal surface could grow ex¬ 
ponentially until disrupt it. 

The situation is more intricate, however, if one keeps 
the perimeter curve fixed, preventing in this way the for¬ 
mulation of an axisymmetric problem from the begin¬ 
ning. By similar arguments, we also expect genus 1 min¬ 
imal area surfaces for small VjLA in this case, but the 
value of the threshold a can be different. Despite of be¬ 
ing unstable as minimal area surfaces, doughnut-shaped 
structures are quite common in many dynamical situa¬ 
tions, ranging from stains left by a coffee droplet [5] to 
extracellular polymeric bacterial coverages [10]. In par¬ 
ticular, toroidal liquid droplets have been obtained by 
exploring some pyroelectric effects on wetting processes 
m- We think our result might be useful to shed more 
light on some of these problems. 


can be expressed by means of elliptical integrals, see, 
for instance, formulas 9 and 7 in the sections 3.152 and 
3.153, respectively, of m- However, the expressions re¬ 
sult rather cumbersome for our manipulations, and we 
chose to solve (A2) numerically. The integrand diverges 
for X = Xmin and x = Xmax, but the divergence is in- 
tegrable and it can be easily circumvented, for instance, 
by introducing the new variable v? = The 

singularity for x = Xmax can be eliminated analogously. 

The equatorial perimeter will correspond to the exter¬ 
nal radius of the doughnut-shaped solution, i.e., to the 
point p* = Ax* such that T'(x*) = 0. The value of x* can 
be determined accurately from (Al) by using a Newton- 
Rapson scheme. The equatorial perimeter will be given 

by 


L = 2tt^. (A3) 

A 

The volume enclosed by the doughnut-shaped surface will 
be given by 


47r r* 

V = ^ J sF{s)ds, (A4) 

and its surface area reads 


Appendix A: Numerics 

It is more convenient for our purposes to introduce 
some dimensionless quantities. From (12), one can intro¬ 
duce the dimensionless function F(x) 


f{p) = jF{Xp), 


(Al) 


^ Stt 




ds. 


(A5) 


The Area x Volume diagram of ^ is constructed from 
(A4) and (A5) by varying d while keeping L given by 
(A3) fixed. 


with 


F{x) = 


d-, 




ds. 


(A2) 


It is clear from (Al) that the constant A can be absorbed 
by a global rescaling of the problem. The integral (A2) 
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